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Abstract 

The kinetics and thermodynamics of structural transformations under pressure de- 
pend strongly on particle size due to the influence of surface free energy [IH4]. By 
suitable design of surface structure |5j, composition [6 J and passivation |2j it is possi- 
ble, in principle, to prepare nanocrystals in structures inaccessible to bulk materials 
[7]. However, few realizations of such extreme size-dependent behavior exist |8j. Here 
we show with molecular dynamics computer simulation that in a model of CdSe/ZnS 
core/shell nanocrystals the core high pressure structure can be made metastable un- 
der ambient conditions by tuning the thickness of the shell. In nanocrystals with 
thick shells, we furthermore observe a wurtzite to NiAs transformation, which does 
not occur in the pure bulk materials. These phenomena are linked to a fundamental 
change in the atomistic transformation mechanism from heterogenous nucleation at 
the surface to homogenous nucleation in the crystal core. Our results suggest a new 
route towards expanding the range of available nanoscale materials. 
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At thermodynamic equilibrium, matter adopts the form that minimizes its total free 
energy [9]. Close to first-order phase transitions, however, metastability of the competing 
phases is often observed; liquid water can be cooled many degrees below its freezing point, 
magnets can withstand oppositely oriented magnetic fields, and diamonds do not transform 
to graphite at ambient conditions. How far one can push a system out of its equilibrium 
phase depends on the microscopic transformation mechanism that determines the height of 
the free energy barrier separating the competing phases. 

Exploiting the metastability of different solid phases is a possible route to creating mate- 
rials with new properties. However, many crystal structures form only when high pressure is 
applied and are unstable under ambient conditions in the bulk. In nanocrystals, on the other 
hand, phase diagrams and microscopic transformation mechanisms can depend strongly on 
particles' size and shape [51 [TO]. The wurtzite to rocksalt transformation in CdSe nanocrys- 
tals, for example, shows an increasing thermodynamic transition pressure and a decreasing 
activation enthalpy with decreasing particles size [21 dH 12] ■ While it is in principle possible 
to extend the metastability of high-pressure structures to ambient conditions by engineering 
the surface properties of nanoparticles, significant insight into the underlying microscopic 
transformation dynamics is required. 

A particularly interesting surface modification is realized in core/shell nanocrystals [T31 . 
where the core material is epitaxially overgrown with a material of identical crystal struc- 
ture [TH - U6] . While the optical qualities of these hetero-materials are well-studied [17] . little 
is known about their structural properties [18]. In modern synthesis methods, materials 
with a lattice mismatch of up to 11% can be combined to form a pristine core/shell in- 
terface [17J. The resulting lattices of both core and shell experience a strong strain that 
depends sensitively on particle size and has the potential of introducing dramatic changes 
to the nanoparticle's structural and kinetic behavior under pressure. 

In this Letter, we report the simulation of spherical wurtzite CdSe nanocrystals of 3 nm 
diameter (510 atoms), that have been epitaxially passivated with ZnS shells of thicknesses 
up to 2.1 nm (5 monolayers). The largest of these core/shell crystals consists of 6918 atoms. 
The particles are modeled with empirical pair potentials designed to reproduce a number of 
properties of the bulk materials [HI [20]. In our simulations, a single crystal is immersed in 
a pressure bath of ideal gas particles [211 E2] at a temperature of 300 K and the pressure 
is increased in steps of 0.2 GPa every 10 ps. These pressurization rates are many orders of 
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magnitude larger than in experiments using diamond anvil cells but are comparable to recent 
shock-wave experiments on CdSe nanocrystals [23] ■ When a pressure of 20 GPa is reached 
after 1 ns, the pressure is released again at the same rate. After reaching ambient pressures, 
the crystals are simulated for another 1.2 ns. We monitor the evolution of the crystal 
structure by calculating atom-coordination numbers based on the radial pair distribution 
functions of core and shell atoms. 

The effect of the ZnS shell on the structure of the CdSe core is dramatic. In Figure [T]B 
we plot the density of the wurtzite core of crystals with different shell sizes as a function of 
external pressure. The density of the core increases significantly with increasing shell thick- 
ness. For a 2 nm shell, this compression effect is equivalent to an additional external pressure 
of 6 GPa, as illustrated in Figure [Tp. This pressure is much higher than the bulk coexistence 
pressure of 2.4 GPa of our CdSe model, and high enough to cause spontaneous transforma- 
tion in bare CdSe crystals. Similarly high pressures were found at the core/shell interface in 
experiments of CdS/ZnS nanocrystals [24J. One might therefore expect the transformation 
in core/shell crystals to happen at lower pressures compared to bare CdSe nanocrystals. 
Quite to the contrary, the upstroke transformation pressure of the core increases strongly 
with increasing shell size, as illustrated in Figure [2] While pure CdSe nanocrystals transform 
at around 6 GPa, transformation pressures of up to 18 GPa are observed for crystals with 
thick shells. 

Our simulations suggest that the increase in upstroke transformation pressure with par- 
ticle size is caused by an increase in thermodynamic transition pressure and a concurrent re- 
moval of favorable surface nucleation sites. We estimated the phase diagram of the nanopar- 
ticles by calculating the pressures at which crystals in different phases have equal enthalpy 
(see Figure [2]). In particular, we consider three combinations of core/shell crystal structures: 
wurtzite/wurtzite, rocksalt/wurtzite, and rocksalt/rocksalt. A shell thickness of 0.5 nm, cor- 
responding to a single monolayer of ZnS, is enough to raise the thermodynamic transition 
pressure by 2 GPa. With increasing shell size, the phase boundary of the all-rocksalt phase 
approaches the bulk thermodynamic transition pressure of our ZnS model. At larger shell 
sizes (3-4 monolayers), the rocksalt/wurtzite phase, featuring distinct crystal structures in 
the core and shell, appears as a stable intermediate between the homogenous phases. The 
thermodynamic transformation pressure from the all-wurtzite to the heterogenous phase is 
fairly insensitive to shell thickness, indicating that a large-shell regime has been reached. 
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The appearance of the rocksalt/wurtzite phase is accompanied by a fundamental change 
of nucleation mechanism. Computer simulations of pure CdSe nanocrystals have shown 
that crystals with low energy surface facets always transform via nucleation events on the 
surface. [U [12] While we observe similar surface nucleation in core/shell crystals with shell 
thicknesses up to 3 monolayers, for crystals with thick shells nucleation happens in the core, 
as illustrated in Figure [3j Despite the nearly constant thermodynamic transition pressure in 
this size-regime, the observed transformation pressures still increased. This indicates that 
the change from heterogeneous to homogeneous nucleation results in an increase in the size 
of the nucleation barrier. 

While a typical core transformation event lasts no longer than 10 ps, transformations of 
shells proceed in steps, creating only confined regions of rocksalt at a time. The hysteresis 
curves in Figure [3] manifest such dynamics. By 20 GPa, however, most shells have completed 
the transformation to rocksalt By comparison, we found that a pure 4 nm ZnS nanocrystal 
remained in the wurtzite structure when subjected to the same pressure protocol, indicating 
that the core also influences the shell. Upon release of pressure, all crystal shells undergo a 
back-transformation: Thick shells transform back to a mixture of wurtzite/zinc-blende, thin 
shells to predominantly amorphous four-coordinated structures. 

Not all crystal cores, on the other hand, underwent a back-transformation. While crystal 
cores with thin shells (< 1 monolayers) and thick shells (> 4 monolayers) transform back to 
mixtures of wurtzite and zinc-blende structures, cores in a broad range of intermediate shell 
sizes remain in the rocksalt structure down to zero pressure, as illustrated in Figure |4j This 
observation suggests a significant increase of the down-stroke kinetic barrier with increasing 
shell thickness. The nature of this increased barrier height is related to the interface between 
the rocksalt core and the re-transformed shell. During the back-transformation of the shell, 
bonds between core and shell atoms stay mostly intact. The re-transformed shell neatly 
passivates the rocksalt core and thus removes favorable nucleation sites at the core surface. 
The amorphization of shells thinner than 1 nm upon back-transformation is a consequence 
of this passivation and the mismatch between the rocksalt and wurtzite structures in the 
core and shell, respectively. Very thick shells, however, transform back into crystalline 
wurtzite/zinc-blende mixtures and induce larger stress at the core/shell interface, facilitating 
the back-transformation of the core. Our simulations demonstrate that rocksalt cores can 
persist on the nanosecond timescale; recent experiments suggest substantial metastability 
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even on experimental timescales [IS] . 

Phase transitions that occur far from equilibrium do not necessarily lead to the phase with 
the lowest free energy. In fact, Ostwald's step rule predicts that a system will transform from 
a metastable phase to the phase with the smallest free energy difference. In an unexpected 
realization of this rule of thumb, we observed a wurtzite to NiAs (B8) transformation in a 
few crystals with thick shells. Figure [5]A. shows three snapshots of a 1.9 nm shell crystal. In 
the course of the transformation a grain boundary between the expected rocksalt structure 
and the NiAs structure builds up in the core and later propagates into the shell at higher 
pressures. Like rocksalt, the NiAs structure is 6-coordinated. Cations are in a rocksalt-type 
coordination environment, while anions are coordinated by a trigonal prism of cations (Fig- 
ure [5^3). The occurrence of the NiAs structure is surprising, since it has not been observed 
experimentally in the pure materials, neither in the bulk nor in nanocrystals. Figure [5p 
shows a plot of the bulk enthalpies of the core and shell materials in the wurtzite, rocksalt 
and NiAs structures as a function of pressure. Throughout the pressure range studied here 
(0-20 GPa), NiAs is never enthalpically most stable. However, it becomes metastable with 
respect to the wurtzite structure at pressures larger than 4.5 GPa and 16 GPa, for CdSe 
and ZnS, respectively. (The transformation illustrated in Figure [5p occurred at 17 GPa.) 
Interestingly, in a recent pressure study of ZnS/CdSe core/shell nanocrystals an unexpected 
Raman peak was observed after the transformation had happened [25] . 

In summary, we have shown that both the kinetics and thermodynamics of the wurtzite 
to rocksalt transformation in CdSe/ZnS core/shell crystals are strongly affected by the 
thickness of the shell. A strong increase in thermodynamic transition pressure with increas- 
ing shell thickness is accompanied by a substantial broadening of the hysteresis, rendering 
the transformed rocksalt cores metastable at ambient conditions. The up-stroke nucleation 
pathway changes from heterogenous nucleation on the surface to homogenous nucleation in 
the core. In thick-shell crystals, the greatly increased up-stroke transformation pressure can 
lead to nucleation of the NiAs structure, which is not observed in the two pure materials. 

The unexpected occurrence of a new high-pressure NiAs structure suggests that other ma- 
terials might be susceptible to a similar phenomenon. By artificially increasing the pressures 
at which solid-solid transformation take place, transformation routes to other, previously un- 
observable crystal structures might become available. Potentially, such an increase can be 
achieved by blocking favorable nucleation pathways through suitable surface modifications, 
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or by using high pressurization rates as obtained in shockwave experiments. 



Methods 

The interatomic pair potential for ZnS was obtained using the same procedure outlined 
in Ref. [20J and consists of a Lennard- Jones part and electrostatic interactions. Potential 
parameters are given in Table [Tj Standard mixing rules were used for interactions between 
unlike atom types. The equations of motion were integrated with the velocity Verlet algo- 
rithm and a time step of 2 fs. The ideal gas barostat was parameterized as in Ref. [22J . 
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17998.4 
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1.18 
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Se 
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TABLE I. Potential parameters for ZnS and CdSe [20] . 

To generate core/shell nanocrystals, we place Cd, Se, Zn, and S atoms on a contiguous 
subset of the lattice sites that define a bulk wurtzite crystal with lattice constants a = 4.3 A 
and c = 7.0 A. A roughly spherical shape is obtained by including in this subset only lattice 
sites lying within a distance r = 15 A + r s h e ii away from the center of a particular (arbitrary) 
unit cell. Surface energies are reduced by subsequently excluding lattice sites possessing 
fewer than two nearest neighbors. Atoms are positioned on the resulting subset of lattice 
sites as follows: (i) Cd and Se are placed on all "a" and "b" sites, respectively, within a 
radius of 15 A of the center, (ii) Zn and S atoms are placed on all remaining "a" and "b" 
sites, respectively, (iii) In a final step designed to mimic effects of CdSe core relaxation prior 
to shell growth, all Cd and Se atoms bound to fewer than two Se and Cd atoms, respectively, 
are replaced by Zn and S atoms, respectively. 

Coarse-grained atom densities were computed as 

p W = (2™y 3/2 Ee*p(^Y (i) 

i=l ^ ' 
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where r« is the position of atom i, and N is the total number of atoms. The effective atomic 
width a = 2 A was chosen to balance smoothness and resolution of p(r). The density 
profile p(r) plotted in Figure [T]A. was computed by averaging p(r) within thin spherical 
shells centered on the nanocrystal's center of mass: The core density values in Figure [T]B 
were further averaged over distances within 8 A of the nanocrystal's center of mass. 

Points of equal enthalpy in Figure [2] were determined in NPT-Monte Carlo simulations 
as described in Ref. [12]. Initial configurations were extracted from the molecular dynamics 
trajectories in which the transformations were observed. Bulk enthalpies in Figure [5] were 
calculated in NPT-Monte Carlo simulations of approximately 500 periodically replicated 
atoms. Ewald sums were used for the calculation of electrostatic interactions [26] . 

Hysteresis curves in Figures [3] and [4] show the fraction of 4- and 6-coordinated atoms in 
the core and shell, respectively. These fractions were calculated as follows: For every con- 
figuration along a given trajectory, the radial pair distribution function g{r) was calculated 
separately for core and shell atoms and averaged over 5 ps intervals. The first minimum of 
g(r) was found numerically and the respective distance used as a cutoff for the calculation 
of the number of nearest neighbors. To minimize interfacial contributions, only core atoms 
within 10 A of the center of mass were included in the calculation; likewise, shell atoms 
within 18 A of the center were excluded. 
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FIG. 1. The ZnS shell compresses the CdSe core. (A) Coarse-grained atom density of a 
1.5 nm shell-crystal at zero pressure as a function of distance r from the center of mass. The different 
densities of the core and shell materials are well visible. (B) Density of the core as a function of 
pressure, for nanocrystals with different shell thickness (legend values indicate shell thickness in A). 
The density of bulk CdSe, obtained from constant pressure Monte Carlo simulations, is shown for 
reference. The dashed line is a fit of the bulk data to the Murnaghan equation of state. Note that 
the sudden increase in density observed for and 0.2 nm shell-crystals is a signature of the wurtzite 
to rocksalt transformation. (C) Core density at zero pressure as a function of shell thickness. The 
right hand ordinate shows the pressure necessary to achieve equivalent densities in bulk CdSe. 
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FIG. 2. Size-dependent transformation pressure. Core upstroke transformation pressures 
(solid black circles) are plotted as a function of shell thickness. At these pressures, the fraction of 
six-coordinated atoms exceeds 0.1 for the first time. The thermodynamic transition pressure of bulk 
CdSe (dashed line) and bulk ZnS (dotted line) are shown for reference. Points of equal enthalpy 
(blue squares and red diamonds), obtained from constant pressure Monte Carlo simulations, give 
an estimate of the nanocrystal phase diagram as a function of shell thickness. The three phases 
are illustrated below the graph as cross sections of a 2 nm shell crystal: wurtzite core and shell 
(WZ), a rocksalt core in a wurtzite shell (RS/WZ), and rocksalt core and shell (RS). 
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FIG. 3. The nucleation mechanism changes with increasing shell thickness. (A) (Left) 
Fraction of 6-coordinated atoms as a function of pressure in the core (black) and shell (red) of 
a 1.2 nm shell crystal. The transformation of both core and shell start around 11 GPa. (Right) 
Cross sections highlight stages of the nucleation process, as seen along the wurtzite c-axis. The 
dashed line marks the interface between core and shell. Atoms that have undergone a change of 
coordination are shown opaque. (For clarity, only clusters of 10 atoms or more are shown.) The 
transformation nucleates at the crystal surface and propagates inward. (B) (Left) Fraction of 6- 
coordinated atoms of a 2 nm shell crystal. The transformations of the core and shell happen at 
different pressures, around 11 and 14 GPa, respectively. (Right) Snapshots show that the nucleus 
is located in the core. 
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FIG. 4. Rocksalt metastability at ambient pressure. (A) Fraction of 6-coordinated atoms 
in the core (black) and shell (red) of a 0.9 nm shell crystal as the pressure is increased to 20 GPa 
and then released again. The rocksalt to wurtzite transformation at around 11 GPa is well visible. 
While the shell undergoes the back-transformation at around 5 GPa, the core remains in the rocksalt 
structure even at zero pressure. (B) Fraction of 4- and 6-coordinated atoms in the core, 1 ns after 
completing the pressure cycle. While crystals with very thin (< 0.4 nm) and thick (> 1.7 nm) 
shells transform back, crystals in a range of intermediate shell thicknesses remain in the rocksalt 
structure. (C) Cross sections of a 0.9 nm shell crystal at different points in the pressure cycle, 
viewed along the wurtzite c-axis. (Left) At an up-stroke pressure of 1.6 GPa, both core and shell 
are in the wurtzite structure. (Middle) At 19.6 GPa the crystal is in the rocksalt structure. (Right) 
A nanosecond after completing the pressure cycle, the rocksalt structure in the core persists. The 
shell has transformed back into a predominantly amorphous 4-coordinated structure. 
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FIG. 5. NiAs structure nucleates at high pressures. (A) Time series of cross-sections of a 
1.9 nm-shell nanocrystal undergoing transforming from wurtzite to NiAs/rocksalt. The crystal is 
viewed along the wurtzite c-axis and the same set of atoms is displayed throughout. (Left) First 
stage of nucleation in the core at 17 GPa. (Center) 4 ps later, no 4-coordinated atoms remain in 
the core and a grain boundary between NiAs (upper left part of the crystal) and rocksalt (lower 
right part) is visible. The shell is visibly strained, but still wurtzite. (Right) At 20 GPa, no four- 
coordinated atoms remain; the NiAs grain-boundary spans the entire crystal. (B) Close-up view 
of a patch of CdSe in the NiAs structure, highlighting the different coordination environments of 
Cd (blue) and Se (red) atoms. (C) Bulk enthalpies per atom as a function of pressure for CdSe 
and ZnS in the wurtzite, rocksalt and NiAs structures. Throughout the pressure range studied (0- 
20 GPa), NiAs is never stable. It is metastable with respect to the wurtzite structure at pressures 
larger than 4.5 GPa and 16 GPa, for CdSe and ZnS, respectively. 
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